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ABSTRACT 

We propose a new explanation for the origin of angular momentum in galaxies and their dark halos, 
in which the halos obtain their spin through the cumulative acquisition of angular momentum from 
satellite accretion. In our model, the build-up of angular momentum is a random walk process associated 
with the mass assembly history of the halo's major progenitor. We assume no correlation between the 
angular momenta of accreted objects. The main role of tidal torques in this approach is to produce the 
random tangential velocities of merging satellites. Using the extended Press-Schechter approximation, we 
calculate the growth of mass, angular momentum, and spin parameter A for many halos. Our random 
walk model reproduces the key features of the angular momentum of halos found in ACDM N-body 
simulations: a log-normal distribution in A with an average of (A) w 0.045 and dispersion ax = 0.56, 
independent of mass and redshift. The evolution of the spin parameter in individual halos in this model is 
quite different from the steady increase with time of angular momentum in the tidal torque picture. We 
find both in N-body simulations and in our random walk model that the value of A changes significantly 
with time for a halo's major progenitor. It typically has a sharp increase due to major mergers, and a 
steady decline during periods of gradual accretion of small satellites. The model predicts that on average 
the A of '-^ 1O^^M0 halos which had major mergers after redshift z = 3 should be substantially larger 
than the A of those which did not. Perhaps surprisingly, this suggests that halos that host later- forming 
elliptical galaxies should rotate faster than halos of spiral galaxies. 

Subject headings: cosmology: theory — cosmology: dark matter — galaxies: evolution — galaxies: 
interaction 



1. INTRODUCTION 

Angular momentum is among the most important quan- 
tities determining the size and shape of galaxies, and yet a 
, detailed understanding of its origins remains a missing in- 
gredient in the theory of galaxy formation. Hoyle (1949) 
was apparently the first astrophysicist who discussed the 
source of galaxy rotation in the framework of a theory 

■ of gravitational instability. He explained galaxy rotation 
] as arising from gravitational coupling with the surround- 

■ ing matter. Alternative theories of the origin of galaxy 

■ rotation due to primordial turbulence and vorticity were 
, also discussed at that time (Wieszacker 1951; Gamow 

1952), but such theories were subsequently ruled out by 
the fact that velocities not arising from gravitation decay 
in an expanding universe, and also because of improving 
constraints on velocities derived from cosmic microwave 
background anisotropics. 

Hierarchical clustering of cold dark matter (CDM) (Blu- 
menthal et al. 1984) is now widely believed to describe the 
origin of structure in the universe, and in this scheme the 
angular momentum of dark matter halos — and eventu- 
ally the rotation of galaxies — is usually thought to be 
produced by gravitational tidal torques in the course of 
the growth of rotation-free initial perturbations. This is a 
well-developed (Peebles 1969; Doroshkevich 1970; White 
1984) and widely accepted theory, and there is little doubt 
that the tidal torques mechanism is connected with the 



generation of angular momentum. To some degree the 
model is supported by N-body cosmological simulations. 
For example, if one considers all the mass which at z = 
ends up inside a given halo, then the angular momentum 
of that mass typically increases linearly with time at early 
stages of the collapse, just as predicted by the tidal torque 
theory (e.g., Barnes & Efstathiou 1987). At later times, 
again in agreement with predictions, the growth in angular 
momentum slows down (Sugerman et al. 2000). 

However, existing models and approximations that try 
to implement the tidal torque scenario do not agree in 
detail with the results of fully nonlinear N-body simula- 
tions. The predicted angular momentum of halos in the 
tidal torque model is typically overestimated by a factor 
of ~ 3 compared to the results of simulations, with a large 
scatter of the same magnitude (Barnes & Efstathiou 1987; 
Sugerman et al. 2000) . This over- prediction factor can be 
lessened if one assumes that the angular momentum stops 
its growth earlier than the turnaround moment (Porciani 
et al. 2002) because of nonlinear effects, but the scatter 
still remains very large. The direction of the spin param- 
eter predicted by the theory also has large errors. Lee & 
Pen (2000) give the average error in the direction of the 
spin of about 57°, which agrees with Porciani et al. (2002) 
who find a mean misalignment of 52° with a scatter ~ 35°. 
The source of the errors is still under debate, but it seems 
that the main effect is due to nonlinear effects which are 
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difficult to take into account in the tidal torque model 
(Porciani et al. 2002). 

We also address another limitation of the tidal torque 
picture: predictions for the spin parameter in this frame- 
work arc for all the matter which at z = is found in a 
given halo. At high redshifts this mass is of course not 
in a single halo. In other words, the tidal torque theory 
does not predict the rotation of any particular progenitor 
of the halo in question. In this paper, we will pursue an 
alternative approach, that of tracing the evolution of the 
angular momentum of the most massive progenitor of the 
present-day halo. This alternative approach is much more 
useful for semi-analytic modeling of galaxy evolution (e.g. 
Somerville & Primack 1999). 

The angular momentum can be c;xprc;ssed in terms of 
the dimensionless spin parameter, which is defined as 

Here J is the angular momentum, E is the total energy, 
and M is the mass of a halo. The value of the spin pa- 
rameter roughly corresponds to the ratio of the angular 
momentum of an object to that needed for rotational sup- 
port (e.g., Padmanabhan 1993). For example, a spin pa- 
rameter of A = 0.05 implies very little systematic rotation 
and negligible rotational support. Typical values of the 
spin parameter of individual halos in simulations are 0.02 
to 0.11 (Barnes & Efstathiou 1987; Rydcn 1988; Warren 
et al. 1992; Steinmetz & Bartelmann 1995; Cole & Lacey 
1996; Gardner 2001). The distribution of spin parameters 
in N-body simulations is well described by the log normal 
distribution: 
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The parameters for the log-normal distribution were found 
to be 0.03 < A < 0.05 and 0.5 < cta < 0.7 for stan- 
dard CDM and various variants (e.g. Warren et al. 1992; 
Gardner 2001). For the ACDM cosmology with matter 
density Qq = 0.3, Q\ = 0.7, h — 0.7, and erg = 1, the 
log-normal parameters were found to be A = 0.042 ± 0.006 
and ax = 0.50 ±0.04 (Bullock et al. 2001b). Note that the 
A parameter of the log-normal distribution (2) is not equal 
to the mean of A; rather, (A) « 1.078A for ax — 0.5 — 0.6. 

In order to study the evolution of the angular momen- 
tum of a dark matter halo it is important to find correla- 
tions of the spin parameter with other parameters of the 
halo and with its environment. Correlation in the direc- 
tions of spins of nearby halos or galaxies is interesting as 
an indicator of the strength of large-scale correlations and 
how much the next infalling satellite "knows" about the 
previous one. Because the tidal torques are due to pertur- 
bations in the gravitational potential, one naively expects 
that there are long range correlations in spins of halos and, 
consequently, in the angular momenta of galaxies. Indeed, 
there are claims that such correlations exist in both obser- 
vational data and in N-body simulations. But the corre- 
lations of galactic rotation axes are measured to be small. 
If L{x) is a unit vector in the direction of the spin axis, 
then Pen, Lee, & Seljak (2000) find that the correlation 
function of directions (i(x) • I/(x -|- r)) is less than 0.05 
for distances ^ 3 /i~^Mpc, and is even smaller at larger 
distances. Even at small separations the effect is small 



and is consistent with no correlations. Results of N-body 
simulations also indicate that spin correlations are weak 
at best (Barnes & Efstathiou 1987). Lee & Pen (2000, 
2001) found the spin-spin correlation function of halos is 
w 0.05. Porciani et al. (2002) found no correlation at 
distances larger than l/i~^Mpc at late stages of evolution 
{z = 0). These very small correlations also suggest that 
material which is accreted by a growing halo has very little 
memory: accretion is mostly random. 

The spin parameter of halos in N-body simulations ap- 
pears to be a very stable statistic and has been shown to 
be independent of most physical parameters. No depen- 
dence has been found on the cosmological model (Barnes 
& Efstathiou 1987; Warren et al. 1992; Gardner 2001; 
Lemson & Kauffmann 1999), on halo environment or on 
halo mass (Lemson & Kauffmann 1999). The only corre- 
lation that has been shown to exist is with the time of the 
last major merger: halos that experienced a recent merger 
have larger spin (Gardner 2001). 

In this paper we investigate the amount of angular mo- 
mentum that is brought in by the random uncorrelated 
accretion of satellites. Because there is no systematic ro- 
tation of the incoming satellites, one might naively expect 
that their contribution to the rotation of the accreting ob- 
ject will be small. But at the same time, each act of accre- 
tion contributes a large orbital angular momentum. Thus 
the halo typically acquires a substantial rotation from each 
large merger, which may later be diluted by many uncor- 
related small mergers. The outcome of the process was 
not obvious in advance, and to our surprise we found that 
the resulting rotation is quite significant. 

We begin, in §2, with an analysis of N-body simulations. 
There are two reasons for this. First, we study the spin 
histories of several N-body halos in order to provide a com- 
parative test for our model predictions. Second, we need 
to characterize the distributions of velocities and angular 
momenta of infalling satellites, because these are essential 
inputs for our model. 

In §3 we present our random walk model for the build- 
up of angular momentum in halos. Results of the random 
walk approximation are presented in §4. We discuss our 
results and the general implications of our model in §5. 

2. SPIN PARAMETERS, HALOS, AND ANGULAR MOMENTA 
OF SATELLITES IN N-BODY SIMULATIONS 

2.1. Numerical simulations 

We use three N-body simulations of a low-density fiat 
ACDM model with the following parameters: Oq = 0.3, 
r^A = 0.7, and h = 0.7. The simulations were done with 
the Adaptive Refinement Tree (ART) code (Kravtsov et 
al. 1997). 

The first simulation, Klypin et al. (2001), was done only 
for three halos, but the resolution was extremely high: the 
mass per particle was 1.2 x lO^/i~^M0 and the (formal) 
force resolution was 100/;^ ^pc. Here we used the power 
spectrum normalization as = 0.9. At redshift z = the 
three halos had virial masses (1.1 — 1.5) x lO^^/i~^M0, 
corresponding to ^ 10^ particles per halo. 

In order to find statistics of velocities and orbital an- 
gular momenta of satellites needed for the random walk 
model, we used two additional simulations, which do not 
have as high resolution as the first one but have many 
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Fig. 1. — Three examples of evolution tracks of galaxy-size ha- 
los in N-body simulations. All halos show fast mass growth at 
high rcdsliifts. At that epoch their spin parameters behaved 
very violently, but subsequently they mostly declined as the 
halo masses grew. 




Fig. 2. — The distribution of directions of the orbital angular 
momenta of accreted satellites relative to the direction of rota- 
tion of the central galaxy-sized halos. From top to bottom each 
row corresponds to halos a, b, and c. Left column is for halos 
at z = 1; the right column is for 2 = 0. The angular momen- 
tum of the central halo corresponds to the northern pole (the 
top point in each panel). There is no clear correlation of the 
distribution of the directions of satellites either with redshift 
or with the direction of rotation of the central halo. 



more halos. One simulation used 256'^ dark matter parti- 
cles within a cubic box of comoving size 60/i~^Mpc with 
mass per particle m,p = 1.1 x IQ^hr^h^^MQ and with a 
(formal) force resolution 1.8/i^^kpc. For this simulation 
we used erg = 1.0. This simulation has 15500 dark matter 
halos at redshift z = with masses from 2 x lO^°/i~^M0 to 
~ lO^^/i~^M0. Another simulation used 512^ dark matter 
particles within a cubic box of comoving size 80h~^yipc. 
The mass per particle was nip = 3.2 x 10^h~^h~^MQ. This 
simulation was run only to z = 3. It had 62,000 halos with 
mass larger than 1.3 x lO^^/i^^Mj;. This simulation was 
used to study statistics of large mergers. 

In our simulations the halos are identified by the 
Bounded-Density-Maxima (BDM) algorithm (Klypin & 
Holtzman 1997). For distinct halos (halos that are not 
inside larger halos) the algorithm defines halos as spher- 
ical objects with average virial overdensity. Halos inside 
larger halos are defined as gravitationally bound lumps of 
dark matter. 

2.2. Evolution of halo spin parameters of major 
progenitors 

Figure 1 shows the evolution of the mass and the 
spin parameter of the major progenitors of our three 
~ lO^^/i~^M0 high resolution halos. The spin parameter 
of the halos clearly changes with time. Even for the most 
"quiet" halo b the spin parameter changed by a factor 1.5 
since redshift 3. For the other two halos the change was 
even larger. Note that there is no steady increase of the 
spin with time; instead, there are large sudden changes in 
A, which clearly correlate with jumps in the masses of the 



halos mostly associated with large mergers. The jumps in 
A are of both signs: most large merger events are associ- 
ated with large spin increases, but halo a had a relatively 
large merger event at z w 0.33, which decreased its spin 
by about a factor of two. There is also a tendency of A 
to decline during periods of gradual mass accumulation, 
which is clearly observed in halos b and c at later stages 
of evolution. 

The halos accrete satellites from all directions. Figure 2 
shows directions of orbital angular momenta of satellites 
which are being accreted by the halos at two redshifts. We 
show satellites that at a given redshift are found in a spher- 
ical shell radius between 0.9 and 1.1 times the virial radius. 
Only satellites that move toward the center are considered. 
The distribution is not totally random: the satellite spins 
are clustered. This clustering of orbital spins comes from 
real-space clustering: there are larger satellites that are 
surrounded by smaller ones. When a group of satellites 
falls in to the central halo, they have almost the same 
angular momentum. 

Besides the clustering, there is no discernible pattern 
in the directions of the satellites. For example, halo b at 
z = has a hole around the northern pole (direction of 
rotation of the central halo). This would indicate that or- 
bits of satellites are far from the plane of rotation of the 
central halo. At the same time, halo c has an excess of 
satellites with the same direction of orbital motion. 

These three halos were chosen to reside in a ~ 10 Mpc 
filament bordering a large void. This is done to roughly 
mimic the environment of the Local Group. The filament 
appeared at very high redshifts. The three halos were 
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Fig. 3. — Probability distribution of the ratio of final to initial 
spin parameter Aj/A^ between subsequent stored time-steps, 
for halos with Mvir> 5 x lO^/i^^M© at 2; = that have un- 
dergone a major merger. Here A; is the spin parameter A' 
at the last stored time-step before the merger, and A^ is the 
spin parameter at the next stored time-step. Panels repre- 
sent: top left, aJl halos; top right, A^ < 0.025; bottom left, 
0.025 < A' < 0.055; bottom right, A- > 0.055. The median 
value of \'f/\'i is printed in the upper left corner of each plot. 
Ifalos with low initial spin parameters typically have large in- 
creases in spin in a major merger. 
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4. — Probability distribution of A'j /Aj for halos with 
5 X 10^^ H^^Mq at 2: = 0. The various panels corre- 
spond to halos with different initial spin parameters and differ- 
ent mass ratios. The rows correspond to increasing initial spin 
values from top to bottom: A' < 0.025, 0.025 < A' < 0.055, 
and X'i > 0.055. The columns correspond to increasing mass 
ratios from left to right: Mf/Mi < 1.1, 1.1 < Mf/Mi < 1.25, 
and Mf/Mi > 1.25, where M/ is the mass of the halo in the 
later time-step and Afi is the mass of that halo's most massive 
progenitor in the earlier time-step. The median value of A^ /X'i 
is printed in the upper left corner of each panel. 



always in the filament, growing by merging with smaller 
halos. A simple naive expectation is to see a lump of satel- 
lites in panels in Figure 2. Its direction would indicate the 
direction of the filament and the rotation of the halo rel- 
ative to the filament. We do not find this: there is no 
obvious large single lump in the diagrams. The problem 
with the naive expectation is that one imagines a thin tube 
with the halo inside it. Even qualitatively this is a wrong 
picture. The filament is much wider (1-2 Mpc) than the 
virial radii of the three halos 200 kpc). When the satel- 
lites eventually fall in to the halo, they come from different 
directions in the filament, not from the main ax;is of the fil- 
ament. Because of the finite thickness of the filament, one 
may even expect to find anti-correlation of orbital spins 
of satellites. Satellites moving parallel to the axis bring 
opposite signs of the angular momentum when they move 
above or below the axis. 

Figure 2 also shows that there is almost no correlation 
with rcdshift for the same halo. Patterns are not repro- 
duced with time. Statistical analysis of the directions 
confirms our visual impression. The measured distribu- 
tion is consistent with an uncorrelated random distribu- 
tion. The data still allow some degree of correlation or 
anti-correlation because of the small number of satellites 
and consequent statistical uncertainties. The total num- 
ber of satellites in our simulations was large, about 300 
per central halo at the end of evolution, yet the number 
of satellites in a shell around the virial radius was much 
smaller, about 40 for each halo. With the total number 



of 120 satellites, statistical uncertainties still allow correla- 
tion or anti-correlation on the level of 10-20%. Altogether, 
the evolution of the spin parameter and the directions of 
accreted satellites indicate that the accumulation of the 
angular momentum happens mostly in a random, uncor- 
related fashion. 

We used the second simulation with a large sample of 
halos to study the correlations of jumps in A and in mass. 
The analyses presented here are based on a new structural 
catalog of halos from this simulation, developed as part 
of the study of the relationship of the structure of halos 
to their merging history in Wechsler (2001) and Wechsler 
et al. (2002), where further details can be found. Briefly, 
all separate bound virialized halos with more than 20 par- 
ticles are cataloged at each of 36 stored time-steps from 
this simulation, and NFW fits (Navarro, Frenk, & White 
1996, 1997) are obtained for all halos with more than 200 
particles. These catalogs were then compared in order to 
determine the merging histories of all halos. The statistics 
presented here are for distinct halos, i.e. those which are 
not subhalos of another halo. 

The simulations indicate that spin parameters of indi- 
vidual halos change very substantially during the growth 
of halos: there are sharp increases in A which correlate 
with periods of fast mass increase, while periods of slow 
mass increase tend to correlate with declines in the spin 
parameter. Thus the features illustrated for three halos in 
Figure 1 are in fact quite general. In Figure 3 we present 
statistics of the change in the spin parameter in major 
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mergers, while Figure 4 includes these statistics for a range 
of changes in the mass of the most massive progenitor. In 
these figures, we actually plot the changes not in A but 
rather in A', defined as (Bullock et al. 2001b; Dekel et 
al. 2001) 



V2M™V;i?™' 
where = GMvir/^Rvir is the circular velocity at the 
virial radius, i?vir, for a halo with virial mass Myir (see 
§3.1 for more details). The spin parameters A' and A are 
approximately equal for typical NFW halos (Bullock et 
al. 2001b), but A' is easier to measure in simulations. 

For Figure 3 we identified all halos with two progenitors 
of mass ratio m/M > 1/3 in the previous stored time-step, 
and investigate the distribution of A'^/A- in these merging 
events. On average, we find that A' increases by about 
25% in major merger events (Figure 3, upper left panel), 
though there is a wide distribution. Not surprisingly, the 
change in A' is more severe for cases in which the initial 
A' was low; in these cases an incoming satellite is likely to 
increase the spin regardless of its impact parameter. Ha- 
los which have low initial spin values can thus experience 
jumps of up to a factor of ten in a single merging event. 
When the initial spin is high, however, merger events can 
serve to decrease the angular momentum if they come in at 
opposite direction to the spin. For very high initial values, 
Ay /A^ decreases slightly on average, occasionally as much 
as a factor of five. 

Major mergers are clearly important since they can 
bring in a large amount of angular momentum at once, 
but minor mergers also play a large role in influencing 
halo spin. In Figure 4, we show the distribution of A^/A^ 
for a variety of initial spin parameters and changes in halo 
mass. We see that the spin of halos with very low initial 
spin values is likely to increase with any incoming material, 
even a small amount. For halos with intermediate initial 
spin values, only major mergers are going to have a large 
effect on their spins, and for halos with very large initial 
spin values, any incoming mass is likely to decrease the 
spin. Thus although major mergers produce the largest 
jumps in A', we find that the entire mass accretion history 
is important for its evolution. 

2.3. Velocities and angular momenta of satellites 

In this section we study the infall velocities of satellites 
into larger halos in the simulations in order to provide 
inputs to the random walk model. We show that the dis- 
tribution of infall velocities can be well parameterized by 
a characteristic Gaussian velocity dispersion and a radial 
anisotropy parameter. In addition, we find that the veloc- 
ity anisotropy changes as a function of the mass ratio of 
satellite to major progenitor, m/M: large satellites tend 
have more radial orbits. Since we would like to include at 
least the basic trend with m/M in our model, we make 
the following simplified division: major mergers satisfy 
m/M > 1/3 and minor mergers satisfy m/M < 1/3. We 
characterize the velocity ellipsoids of each of these popula- 
tions separately. We sorted all halos by mass and selected 
pairs of halos starting with the most massive ones. Be- 
cause we are interested in studying satellites at the time 
of accretion, we selected only pairs separated by approxi- 
mately the virial radius Ryir of the larger halo in each pair. 



Three criteria were used: (1) the distance between a pair 
of halos was required to be (1 ± 0.2)i?vir, (2) the relative 
velocity of halos had to be less than 2Vc, where Vc is the 
maximum circular velocity of the larger halo, and (3) the 
halos had to have negative relative radial velocities (cor- 
responding to inward motion), i.e. the angle between the 
line from the smaller to the larger halo and their relative 
velocity vector should be between 0" and 90°. 

The large sample of halos in our 60/i~^Mpc simula- 
tion allows us to study the statistics of relative velocities 
and angular momenta of satellites falling onto larger ha- 
los. We did so by averaging results for expansion factors 
a = 0.972,0.982, and 1.00. Numerical effects — force and 
mass resolution and specifics of our halo finder — result in 
a rather complicated selection of mergers of different mass 
ratios. The limited force and mass resolution result in de- 
struction of most of satellites of small halos. This is the 
so called ovcrmerging problem. Because of this problem 
the galaxy-size halos in this simulation can not be used for 
studies of minor mergers, since their satellites were not re- 
solved. Small subhalos were resolved for large groups and 
clusters, which had dozens or hundreds of them (Colin et 
al. 1999). Thus, we studied minor mergers for "host" 
halos with masses larger than lO^^/i~^M0. Overall, there 
were 543 satellites at radius (1 ± 0.2) of the virial radius of 
the "host" of which 301 were falling onto the "host" and 
the rest were moving out. Merging with large satellites 
was studied mostly at high redshifts. This was done using 
the 60/i~^Mpc simulation at z = 1 and z = 3. We also 
used the 80/i~^Mpc simulation at 2; = 3 to probe statistics 
of satellites with different masses. 

Figure 5 shows several statistics for the minor mergers 
in our catalog. The upper right panel shows the distribu- 
tion of infall angles, ^, of the smaller halos, where <j) is the 
angle between the velocity and radius vectors of the satel- 
lite. The distribution of relative infall velocities in units 
of Vc of the larger halo is shown in the lower left, and 
the distribution of the relative specific angular momenta 
of satellites is shown in the bottom right, plotted in units 
of Lvir = VcRvii of the larger halo, where L = |r x V|. 
The error bars represent Poisson errors. 

In the top left panel we show the distribution of the rel- 
ative velocity as a function of (p. There is a small bias in 
the distribution of velocity vs. angle, with higher infalling 
velocities corresponding to more radial infall directions. 
The distribution of the angles in the top right panel shows 
that there are more infalling halos with radial than with 
tangential motions. 

In general, the magnitudes of the infall velocities are dis- 
tributed around the circular velocity of the massive pro- 
genitor (bottom left panel). The distribution of the veloc- 
ity can be fit with a 3D Maxwell-Boltzmann distribution 
{a, = 0.62Vc, = 3^2): 

Pi"^) = \f^K^M-VV2al). (4) 
V TT crj; 

This fit provides a rough characterization of the distri- 
bution, but does not fully capture its nature. This is 
likely because the the velocities are anisotropic, and the 
Maxwell-Boltzmann distribution describes a system with 
isotropic velocities. The distribution of satellite angular 
momenta (bottom right panel), or, equivalently, tangential 
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Fig. 5. — The distributions of the angles, velocities, and or- 
bital angular momenta of infalling satellites for minor mergers 
identified in the simulation. The angle if) represents the angle 
between the radius and velocity vectors. The top right panel 
shows the distribution of (j), and the top left panel shows the 
distribution with respect to of the infall velocity in units of 
the virial circular velocity of the larger halo. The bottom pan- 
els show the distributions of the relative velocities (left) and of 
the angular momenta (right) . The error bars show the Poisson 
errors. Smooth curves are the fits of the distributions (see text 
for details). 
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Fig. 6. — The dependence of the specific angular momen- 
tum L of a satellite on the ratio of the satellite maximum cir- 
cular velocity to the maximum circular velocity of the host 
halo. The specific angular momentum is measured in units of 
Z/vir = i?virVvir. The spccific angular momentum of small satel- 
lites is larger than that of large ones; the line shows the fit of eq. 
(7). There is no dependence on the redshift. Different mark- 
ers show results for different simulations and redshifts. Stars, 
open circles, and triangles are for the 60/i~^Mpc simulation at 
redshifts z =Q, 1, and 3 respectively. Crosses are for the three 
high resolution halos at z = 0. Solid circles are for z = 3 halos 
in the 80/i~^Mpc simulation. 



velocities, is described reasonably well by a (2D) Maxwell- 
Boltzniann shape: 

p{L) - ^exp(-LV2^i), (5) 

with = 0.59Lvir, {L"^) = 2a|. A good fit is obtained 
because the tangential velocities are fairly isotropic. 

These results suggest that the velocities of satellites 
are well-described by a 3D Gaussian distribution with a 
substantial radial velocity anisotropy. We further have 
checked the distribution of each component of velocity and 
found them to be roughly Gaussian. Based on this, it is 
justified to summarize the velocity distribution of infalling 
satellites with two parameters defining a velocity ellipsoid, 
one describing the velocity anisotropy, 

and the other describing the radial Gaussian velocity dis- 
persion, ar- For example, in the case of isotropic orbits 
(j\ = 2(7^ and /3 = 0, and in the case of nearly radial 
orbits «C 2cr^ and /3 « 1. 

From the 60 Mpc numerical simulation, the velocity 
anisotropy and tangential velocity dispersion were found 
to depend on the ratio m/M of the mass of the satellite 
to that of the host galaxy. For major mergers (3 = 0.8 
and fT^ ~ 0.45Vc, and for minor mergers /3 = 0.6 and 
i7_L = 0.71Vc. These values are in good agreement with 
Colin et al. (1999). Note that major mergers are sig- 



nificantly more radial than minor mergers and bring in a 
factor of 1.6 less specific angular momentum. 

We also studied the velocities of satellites in the sim- 
ulation of the three galaxy-size halos. The statistics of 
satellites in each halo are much better in this simulation. 
We can also easily track the evolution of each halo with 
redshift. The results are consistent with what we found 
for the larger simulation. For minor mergers the average 
velocity anisotropy was (3 = 0.67 at z = and p = 0.53 
at z = 1, however, the difference in this value is not sta- 
tistically significant, and on average the value remains un- 
changed (/3 = 0.6). The three halos in the first simula- 
tion have masses of ~ lO^^/i~^M0. This is significantly 
smaller than 10^"^ — lO^'^h~^MQ in the case of the large- 
box second simulation where, as we explained above, most 
of the " satellites" were coming from group-size central ha- 
los. Thus, the results for the three halos suggest that there 
is no (significant) trend of (3 with the mass of the central 
halo. 

The dependence of the angular momentum of infalling 
satellites on the mass of the satellites is very important 
for accurate estimates of the spin parameter A. We use 
all available simulations to measure this dependence. Be- 
cause the mass of a satellite is not always a well defined 
quantity (e.g., it depends on an uncertain truncation ra- 
dius and on details of removal of unbound particles), we 
prefer to use a more reliable characteristic — the maxi- 
mum circular velocity. Figure 6 shows the dependence of 
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the virial concentration c = Rvir/rs, where the virial ra- 
dius i?vir is defined by the virial mass Mvir and by the 
parameters of the cosmological modeL 



1/3 



(8) 



(9) 



with (Bryan & Norman 1998) 

Sth ^ (IStt^ + 82x - 39x^)/{l + x), 
X = -(1 - no)ay[no + (1 - f^o)a^], 
and a = (1 + In the ACDM model Sth 



a <C 1 and Sth — 340 at a = 1 
truncated at the virial radius is 

1 



(10) 

(11) 
180 at 

The total energy of a halo 



W = 



cgic) 



Here 



and 



1 CM^ 



/(c)=ln(l+c) 



1 + c 



g{c) = 1-2 



ln(l + c) 
1 + c 



{1 + cf 



The spin is then 



A = 



J 



[C9{c)] 



1/2 



(12) 



(13) 



(14) 



(15) 



the specific angular momentum of a satellite L on the ra- 
tio of the satellite maximum circular velocity ysatcintc to 
the maximum circular velocity of the host halo Vhost- Dif- 
ferent markers show results for different simulations and 
redshifts. Error bars show la shot noise estimates. 

Figure 6 indicates that the specific angular momentum 
of small satellites is larger than that of the large ones. It 
also shows that there is no dependence of L/Lsa.% on the 
redshift. The same type of analysis shows that the radial 
velocity dispersion does not depend on the mass of the 
satellite: (Jr/Vc = 0.7 ±0.1. Thus, from eq. 6, velocity 
anisotropy (3 varies from /? = 0.5 for the small satellites to 
/3 = 0.8 for the major mergers. 

We use the following approximation for the rms of the 
specific angular momentum: 

_ r 0.7, if T4at/ Vhost < 0.4, 

-^Z-^™ ~ \ 0.9-0.5(-^), otherwise. 

Figure 6 shows that this is a reasonable summary of the 
simulation data. 

3. RANDOM WALK APPROXIMATION 

The main idea of our random walk model is that halo 

angular momentum is built up from the summed contribu- 
tions of uncorrelated orbital angular momenta of accreted 
objects during the process of halo assembly. The model 
needs several ingredients. First, we need to have the mass 
accretion history for each halo (a merger tree) that tells 
us how many satellites of what mass are accreted by the 
major progenitor as a hmction of redshift. For this, we use 
the extended Press-Schechtcr (EPS) approximation (Press 
feSchechter 1974; Bower 1991; Bond et al. 1991;Lacey& 
Cole 1993) and the EPS merger-tree method of Somerville 
& Kolatt (1999). ^ Second, we need to know the position 
of each satellite at the moment of accretion. The position 
is assumed to be randomly distributed on a sphere whose 
radius is the virial radius of the main progenitor, with the 
main progenitor in the center. 

Finally, we need to know the velocities of accreted satel- 
lites. The three components of the velocities of each satel- 
lite are assumed to be Gaussian, with a velocity anisotropy 
described by eq. (6). We use cr^ = 0.7T4 regardless of the 
mass of the satellite. The tangential velocity dispersion is 
assigned in the following way. For a given satellite mass 
we estimate its concentration parameter (see eqs. (16-17) 
below). The virial mass and the concentration define the 
maximum circular velocity Ksat- We then use eqs. (6- 
7) to find the tangential velocity dispersion. A realiza- 
tion drawn from the Gaussian distributions gives the three 
components of the velocity of the satellite. 

3.1. Mass, concentration, and spin parameter 

The model described above provides all of the informa- 
tion needed to calculate the spins of an ensemble of halos. 
The spin parameter A is determined by calculating the en- 
ergy for an NFW halo in virial equilibrium ^ and then using 
eq. (1). For an NFW halo, p{r) oc [{r / r s){l + r / r sf]~^ , 

® This method describes fairly accurately the mass accretion history of major progenitors in simulations (Somerville et al. 2000; Wechsler et 
al. 2002). 

Our assumption of virial equilibrium will likely break down just after a major merger occurs. We assume that equilibrium will be restored 
over a short enough timescale that the general nature of spin evolution will not depend sensitively on this detail. 
* We also neglect the correlation between concentration and halo merging history (Wechsler 2001; Wechsler et al. 2002). 



2KA^viri?vir /(c) 

For the median concentration c, we use the results of Bul- 
lock et al. (2001a): 

c(A/f™,z) = a(M™)/(l + z), (16) 

a(Afvir) ^ 54.8 - 3.341og(M™/MQ). (17) 

In this paper we have assumed for simplicity that the 
halo mass and redshift uniquely determines its concentra- 
tion. N-body simulations indicate that there is a scatter in 
c at fixed mass of the order of 35% (Bullock et al. 2001a; 
Wechsler 2001). In principle, we could include those de- 
viations in our random walk model, but we decided to 
neglect them because the deviations in A as a result of 
this scatter are less than 10 percent. ® 

Thus, knowing the mass and the angular momentum of 
a halo one can find the halo's concentration and its spin 
parameter. 

3.2. Random walk model for the angular momentum 

Our main results are derived from two sets of Monte- 
Carlo merger trees for two different final masses, Mmax- 
Each set of runs traced the merging history of the most 
massive progenitor back to 2 = 10 or until its mass 
fell below Mmm- The first set had M^ax = IO^^Mq, 
-^min = 1O^M0, and a mass resolution (minimum mass 
of a satellite) of M^^s = 10^ Mq. The second set of trees 
was done with M^ax = 10^^ Mq, Mmin = 3 x IO^Mq, and 
Mres = 1O^M0 . We used 300 mass tracks for each set and 
each track was run 100 times with different realizations of 
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Fig. 7. — The example of two typical evolution tracks (solid 
and dashed lines) of the two dark matter halos. On the left 
panels the mass (top) and the specific angular momentum 
(bottom) are plotted as functions of the expansion parameter 
a= (1 + z)^^. The right top panel shows the mass ratio of the 
most massive merging satellite to the main progenitor, where 
the case of major merger is defined as mass ratio m/M > 1/3. 
The spin parameter as a function of the expansion parameter 
is plotted in the right bottom panel. 



Fig. 8. — Comparison of the distributions of the spin parame- 
ter at redshift z = 0. The solid line represents the distribution 
of A derived using the random walk model. The distribution of 
A from our 60 h^^ ACDM simulation is shown by the dashed 
curve. The log-normal fit to the random walk model is pre- 
sented by the dot-dashed curve. The parameters of the log- 
normal distributions axe in Table 1. 



satellite velocities. Unless otherwise stated, we use these 
two sets of merger trees. 

The evolution of the mass of two typical 10^'^ Mq tracks 
as a function of the expansion parameter is shown on the 
top left panel in Figure 7. The right top panel shows the 
ratio of the mass of the biggest satellite to the mass of the 
main progenitor. The mass of the main progenitor does 
not change smoothly, because we assume that all new in- 
coming mass instantly merges with the main progenitor. 

When the mass of the merging object is significant, there 
are often dramatic changes in the specific angular momen- 
tum (bottom left panel) as well as in the spin parameter 
(bottom right panel). This also was seen for simulated ha- 
los in §2. Of course, major mergers do not always increase 
the spin parameter; depending on the velocity orientation, 
such mergers can also decrease the halo's spin. In spite 
of all these changes in the spins of individual halos, espe- 
cially in the past when the number of major mergers was 
significant, we will show below that the distribution of the 
spin parameter doesn't depend on time or on the mass of 
the main progenitor. 

These examples show that the angular momentum 
doesn't increase linearly. Whether there are increases or 
decreases in any given merging event depends on the mass 
of the satellite, its direction of motion, and its velocity. 
One can also sec that at the later stages of the evolution 
of a dark matter halo, when the mass of the main progen- 
itor doesn't change much, the angular momentum doesn't 
have dramatic changes either. This behavior is very simi- 
lar to that found in the N-body simulations, although the 
comparison cannot be made precise because of the finite 



differences (up to ^ 0.5 Gyr) between the stored time-steps 
in the simulation compared to the instantaneous approxi- 
mation used in our random walk model. 

4. RESULTS 

4.1. Comparison with numerical simulations 

We begin by illustrating that our random walk approxi- 
mation naturally produces a log-normal distribution in A. 

Figure 8 shows a comparison between N-body results and 
our random walk model of spin acquisition. The simula- 
tion corresponds to ACDM {^q = 0.3, h = 0.7, cth = !)• 
For the random walk model we used 200 Monte Carlo 
merger trees with a final {z = 0) mass of the main pro- 
genitor M = 1O"'^'''M0. Each mass track was run 100 times 
with different realizations of satellite velocities. The distri- 
bution of the spin parameter from the N-body simulation 
is shown by dashed line while the distribution from the 
random walk model is shown by the solid curve. As one 
can sec, there is a reasonably good agreement between the 
N-body results and the random walk model. To demon- 
strate that these distributions are well described by the 
log-normal function we also show a log-normal fit to the 
random walk data (the dashed-dotted curve). Parameters 
of the fit are given in Table 1 , along with spin parameters 
measured from N-body simulations. 

4.2. Dependence of the spin parameter on mass and 

redshift 

As we discussed in §1, N-body simulations indicate that 
the A distribution does not vary as a function of redshift or 
halo mass. Here we test our model against these results. 
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X 

Fig. 9. — The distribution of the spin parameter for halos 
at z—0 and z—3, for halos whose main progenitor has mass 
W^^Mq or IO^'^Mq at redshift z = 0. The spin distribution 
has rather weak dependence on mass and redshift. 



We analyzed two different sets of merger trees: one with 
major progenitor mass M = IO^^Mq at redshift z = 0, 
and the other with mass M = 10^^ Mq (sec §3). For both 
cases the distributions of the spin parameters were exam- 
ined at two different redshifts: z = and z = 3. These 
arc shown in Figure 9. The spin parameter distributions 
derived from the random walk model show no correlation 
with either the redshift or the halo mass; the distributions 
overlap and each has mean value (A) = 0.045. 

4.3. Major mergers and the spin parameter 

As was just shown, the spin parameter distribution 
doesn't depend on the redshift or on the mass. But Figure 

7 indicates that A of individual halos in the random walk 
model is sensitive to whether they have recently had major 
mergers, as we found in N-body simulations (Figures 1,3). 
In this subsection we study the influence of these types of 
mergers on the distributions of the spin parameter. 

Wc study the effect of major mergers on the distribu- 
tion of the spin parameter by splitting all halos into two 
groups: the halos that didn't have any major mergers in 
the past roughly 11 billion years (0 < 2; < 3), and the 
halos that did have a major merger during this period. 
Since major mergers destroy galactic disks and produce 
spheroidal stellar systems (see e.g. Barnes 1999), while 
mergers with mass ratio less than 1/10 probably do not 
(e.g. Walker, Mihos, & Hernquist 1996), and stellar disks 
are typically up to about 10 Gyr old, only the halos that 
did not have a major merger since z ~ 3 could host mature 
spiral galaxies. The distributions of A for these two groups 
of halos are shown in Figure 10. The mean value of the 




0.05 0.1 0.15 0.2 



X 

Fig. 10. — The distribution of the spin parameter for dark 
matter halos that had (dashed line) and didn't have (solid line) 
a major merger event since redshift z = 3. Final mass of the 
main progenitor at = is IO^^Mq. The halos that had major 
mergers during the last w 11 Gyrs on average have larger (by 
25%) spin parameter than those which did not experience a 
major merger. 



spin parameter is 0.0465 and 0.0569 for no major merg- 
ers and major mergers since z = 3 respectively, in other 
words, recent major mergers increase the spin parameter 
almost by 25%. Thus, we find the same result that was 
found in the N-body simulations: the spin parameters are 
sensitive to major mergers. 

5. DISCUSSION 

We have presented a new model for the origin of the 
angular momentum of dark matter halos in which angular 
momentum is built up in a random walk fashion by mass 
accretion events. The evolution of halo angular momen- 
tum in this picture is quite different from that inferred 
from the standard tidal torques argument, in which the 
angular momentum grows steadily at early times and its 
growth flattens out at late times. Unlike this standard pic- 
ture, we do not track the angular momentum of all matter 
that is in the present-day halo; we track only the major 
progenitor of the halo. The evolution of the spin parame- 
ter A of the major progenitor is a random process in which 
the the value of A can vary by factors of a few over its 
history, with a tendency for sharp increases due to ma- 
jor mergers and steady decline during periods of gradual 
accretion of small satellites. These general properties of 
evolution, including large fluctuations in A, are confirmed 
by N-body simulations. 

Our model predicts that the distribution of the spin pa- 
rameter is well approximated by the log-normal distribu- 
tion, with mean and dispersion that do not depend on red- 
shift or on the mass of the main progenitor. The model 
predicts that on average the spin parameters of halos that 
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Table 1 

Parameters of log-normal approximations for the spin distribution in a ACDM Model. 



Model 




A 


<J\ 


Comments 


Numerical simulations: 










all halos at z = 




0.0445 


0.563 


Figure 8 


halos that had no major merger since z - 


= 3 


0.034 






halos that had a major merger since z = 
Random walk model: 


3 


0.044 






halos with M = IO^^Mq, z = 




0.0465 


0.556 


Figure 8 


halos with M = 10^^ Mq & 10^^ Mq, z = 


& 2 = 3 


0.0450 


0.560 


Figure 9 (mean values) 


halos that had no major merger since z = 


= 3 


0.0465 


0.333 


Figure 10 


halos that had a major merger since z = 


3 


0.0569 


0.407 


Figure 10 



had major mergers after redshift z = 3 should be con- 
siderably larger than the spin parameters of those halos 
that did not. Perhaps surprisingly, this implies that halos 
which host later-forming elliptical galaxies should rotate 
faster than halos of mature spiral galaxies. 

There may actually be observational evidence of this. 
Planetary nebulae in the two elliptical galaxies NGC 5128 
(radio continuum source Cen A) (Hui et al. 1995) and 
NGC 1316 (Fornax A) (Arnaboldi et al. 1998) are rotating 
rapidly about their minor axes at large radii. The metal- 
rich globular clusters of NGC 5128 show rotation similar to 
that of the planetary nebulae. NGC 5128 is regarded as a 
prototypical remnant of the merger of two disks, and NGC 
1316 is also considered a likely remnant of a fairly recent 
merger. From numerical simulations (Barnes 1992; Hern- 
quist 1993) of mergers between disk galaxies with bulges, 
one indeed expects the angular momentum of the remnant 
to be concentrated in the outer regions. Similar studies of 
the planetary nebulae of other elliptical galaxies such as 
NGC 4697 detect less evidence of strong rotation (Mendez 
et al. 2001), and stellar spectra of only a relatively small 
fraction of elliptical galaxies show strong rotation at large 
radii (with Vr/(T > 1/2), e.g. NGC 1395 and NGC 1604 
(Franx, lUingworth & Heckman 1989). The globular clus- 
ters of M87 and a few other giant elliptical galaxies are 
also rotating rapidly at large radii (Kissler-Patig & Ger- 
hardt 1998; Cohen 2000; Cote et al. 2001). Galaxies 
like the cD M87 are thought to form from multiple merg- 
ers, and simulations of this process predict that while ve- 
locity dispersion should dominate in the central regions, 
rotation can grow to ?;,./cr ^ 1 in outer regions (Weil & 
Hernquist 1996). Radial velocity measurements for the 
planetary nebulae and globular cluster systems of nearby 
elliptical galaxies should become much more feasible with 
the advent of multi-object spectrographs on 8-meter class 
telescopes. Once data are available on a larger number 
of such systems, it will be interesting to see whether the 
statistics are consistent with predictions of the random 
walk model for the origin of galaxy rotation that we have 
proposed here. 

Implementation of the random walk model is a rather 
complicated procedure if one wants to reproduce N-body 
results with an accuracy better than about 20 percent. It 
requires measurements of many parameters of halos and 
satellites (e.g., halo concentrations for different redshifts 
and merging trees). The most important and, unfortu- 



nately, the most uncertain are the parameters of the ve- 
locity ellipsoid of accreted satellites for different masses of 
the satellites. In principle, these parameters could vary as 
a function of redshift as well. Estimates for the numerous 
minor mergers seem to be reasonably reliable (sec §2.3). 
Major mergers are rarer and therefore the estimates are 
more uncertain. Generally we find that the trajectories of 
mergers with high m./A'I ratios are more radial than their 
low m/M counterparts, but it is difficult to measure the 
parameters of the velocity ellipsoid accurately. 

The model assumes that there is no net angular mo- 
mentum of infalling satellites and that the satellites are 
accreted in random fashion. This picture is supported by 
results of N-body simulations (§2.2). Yet the results allow 
some degree of either correlation or anti-correlation of or- 
bital spins of the satellites. The allowed effects are rather 
small (less than 10-20%), but they still may be important 
for accurate predictions of the angular momenta of halos. 

Our colleagues Mailer, Dekel, & Somerville (2002) have 
investigated a simplified version of what we call the ran- 
dom walk model for the origin of the angular momentum 
of dark matter halos, in comparison with a version of the 
tidal torques model more sophisticated than the standard 
one discussed briefly in §1 of the present paper. They 
find that both models can reproduce the log-normal dis- 
tribution of halo spin parameters seen in simulations, with 
appropriate tuning of the model parameters. Their work 
is thus complementary to that presented here. 

We expect that our random walk model for the origin 
of halo angular momentum through the accretion of satel- 
lites will be useful both conceptually and practically. In 
addition to giving a much better fit to the distribution 
of halo spin parameters than the tidal torque theory, it 
also accords better with the hierarchical nature of the pro- 
cess by which halos grow in the cold dark matter model. 
Moreover, because it is based on the same extended Press- 
Schechter formalism used in semi-analytic approaches and 
because it is so simple to implement, we expect that it 
will prove invaluable in improving the treatment of angu- 
lar momentum in modeling the formation and evolution of 
galaxies. For these same reasons, it will be worthwhile to 
test and improve the random walk model. 

We conclude with some comments about the implica- 
tions of our model for disk galaxy formation and the an- 
gular momentum problem in disks. This random walk pic- 
ture is somewhat at odds with the formalism made popular 
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by Fall & Efstathiou (1980) in which the gas is initially well 
mixed within a smoothly rotating halo, and subsequently 
falls in to form an angular momentum supported disk.^ In 
our random walk model, the halo angular momentum ob- 
tains its smoothly varying distribution through a scries of 
clumpy merger events and subsequent relaxation. Unlike 
a picture in which tidal torques determine the spin of a 
halo at very early times, when the gas and dark matter 
was well-mixed, our scenario would suggest that the an- 
gular momentum distribution of gas could be considerably 
different than that of the dark matter. The merging satel- 
lite galaxies are likely to be strongly affected by the tidal 
field of the central object, and their gas is likely to shock 
against the gas in and around the host. Since it is known 
that the angular momentum distribution in dark halos is 
quite unlike that in disk galaxies (Bullock et al. 2001b; 
van den Bosch, Burkert, & Swaters 2001), generally with 
an excess of low-spin material, our random walk picture 



might provide an interesting avenue for exploring solutions 
of this angular momentum distribution problem and also 
of the well-known angular momentum overcooling problem 
(e.g. Navarro & Steinmetz 1997). Work in this direction 
is under way. 
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